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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted,  to  metric  (Si)  units  as  follows: 

Multiply Bv  To  Obtain 


feet  0. 
miles  (U.  S.  statute)  1. 
square  miles  (U.  S.  statute)  2. 
cubic  feet  0. 
cubic  feet  per  second  0. 
feet  per  second  0. 


301*8 

metres 

60931+1+ 

kilometres 

589988 

square  kilometres 

02831685 

cubic  metres 

02831685 

cubic  metres  per  second 

301+8 

metres  per  second 

NUMERICAL  ANALYSIS  OF  TIDAL  CIRCULATION  FOR  LONG  BEACH  HARBOR 


EXISTING  CONDITIONS  AND  ALTERNATE  PLANS  FOR  PIER  J 
COMPLETION  AND  TANKER  TERMINAL  STUDY 

PART  I:  INTRODUCTION 

Objective 

1.  The  purposes  of  this  phase  of  the  Long  Beach  Harbor  model 
study  were  to  numerically  investigate  tidal  circulation  in  existing 
basins  and  to  define  and  evaluate  the  impact  that  alternate  plans  for 
Pier  J completion  for  oil  tanker  berthing  and  general  cargo  facilities 
would  have  on  existing  harbor  circulation. 

S~.  This  report  describes  the  verification  of  the  numerical  model 
using  prototype  data  and  physical  model  data  for  existing  conditions, 
the  evaluation  of  the  model's  predictive  capability  by  applying  it  to  a 
harbor  configuration  which  had  been  previously  studied  in  a physical 
model  of  Los  Angeles  and  Long  Beach  Harbor,  and  the  use  of  the  model 
in  a completely  predictive  mode  to  evaluate  the  effect  of  alternate 
Pier  J configurations  on  tidal  circulation.^-^ 

Background 

3.  Los  Angeles  and  Long  Beach  Harbors  are  adjacent  ports  in  San 
Pedro  Bay,  California.  They  constitute  separate  political  entities  and 
are  administered  by  separate  port  authorities.  Modifications  to  the 
existing  harbors  have  been  proposed  by  both  port  authorities,  the 
U.  S.  Army  Corps  of  Engineers,  and  the  city  of  Long  Beach.  Proposed 
modifications  are  so  extensive  that  a careful  examination  of  their 
effects  on  the  existing  harbor  and  expanded  harbor  facilities  is  re- 
quired to  reduce  the  possibility  of  undesirable  effects  that  could  prove 
irreversible  or  expensive  to  correct.  Therefore,  Congress  directed 
the  Corps  of  Engineers  to  build  a physical  model  of  Los  Angeles  and 


Long  Beach  Harbors  and  to  conduct  studies  of  the  harbors  and  proposed 
modifications.  Construction  of  the  model  at  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  began  in  January  1973  and  was  com- 
pleted in  August  of  that  year. 

4.  An  extensive  prototype  data  collection  program  (described  in 
detail  in  Reference  l)  was  performed  to  provide  data  for  verification 
of  the  physical  model  and  to  describe  existing  conditions  in  the 
harbors.  The  physical  model  was  verified  using  prototype  data,  and 
base  circulation  tests  were  conducted  as  reported  in  Reference  2. 

In  addition,  studies  of  several  alternate  harbor  configurations  have 
been  conducted  in  the  physical  model.  Of  particular  interest  relative 
to  the  present  study  is  the  harbor  configuration  referred  to  as 
plan  1A-2. ^ 

5.  In  November  1975,  the  Port  of  Long  Beach  requested  WES  to 
conduct  a study  of  the  effects  on  tidal-generated  circulation  in  Los 
Angeles  and  Long  Beach  Harbors  produced  by  proposed  modifications  of  the 
Pier  J complex.  Primarily,  because  of  time  constraints  and  scheduling 
problems  associated  with  the  physical  model,  it  was  decided  to  conduct 
the  study  using  a two-dimensional,  depth- aver aged,  numerical  model  of 
the  hydrodynamic  equations.  The  numerical  model  would  be  verified  using 
physical  model  data  for  (a)  the  present  harbor  configuration  and  (b) 
for  a proposed  modification  to  the  harbor  known  as  plan  1A-2.  After 
model  verification,  the  numerical  model  would  be  used  to  predict  the 
effects  on  tidal-generated  circulation  of  proposed  modifications  to  the 
Pier  J complex.  A Pier  J modification  is  desired  that  will  satisfy  the 
requirements  of  a tanker  terminal  and  at  the  same  time  produce  a mini- 
mal overall  circulation  change  from  those  conditions  existing  in  the 
present  harbor  configuration. 
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PART  II:  THE  NUMERICAL  MODEL 


Formulation  of  the  Model 


6.  A complete  mathematical  description  of  the  hydrodynamic  flow 
in  a harbor  or  estuary  would  require  that  the  velocity  and  density  be 
completely  specified  for  every  point  in  the  system  at  all  times: 

u = u(x,y,z,t)* 

p = p(x,y,z,t) 


where 

x = longitudinal  coordinate  measured  along  the  estuary  axis 
y = transverse  coordinate 
z = vertical  coordinate 
t = time 

Unfortunately,  with  existing  computers,  a generalized  three-dimensional 
model  with  the  needed  time  and  spatial  resolution  is  not  available. 
Because  of  the  difficulties  in  formulation  and  executing  a three- 
dimensional  program,  researchers  have  devised  a variety  of  numerical 
models  of  various  degrees  of  simplification.  In  general,  the  simpler 
the  model  the  less  reliable  and  less  adaptable  it  is  to  changing 
conditions. 

7.  A two-dimensional  approach  that  produces  a pseudothree- 
dimensional effect  was  utilized  in  the  numerical  investigation.  The 
vertical  components  of  velocity  and  acceleration  are  neglected  and 
the  general  three-dimensional  governing  hydrodynamic  equations  are 
integrated  over  the  water  depth.  A depth-averaged  two-dimensional 
flow  field  is  obtained  but  three-dimensional  geometry  can  be  con- 
sidered. The  most  important  approximations  used  in  the  model  are 


* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  A). 
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where 

u = depth-averaged  velocity  component  in  the  x direction 
t = time 

x,y  = rectangular  coordinate  variables 

v = depth-averaged  velocity  component  in  the  y direction 
g = acceleration  due  to  gravity 

H = water  level  displacement  with  respect  to  datum  elevation 
f = Coriolis  parameter 

Rx,R  = the  effect  of  bottom  roughness  in  c and  y directions 

L ,L  = the  acceleration  effect  of  the  wind  stress  acting  on 

X V . 

J the  water  surface  m the  x and  y directions 
h = water  depth 

9.  The  continuity  equation  has  been  obtained  by  integrating 
across  the  water  depth  and  applying  kinematic  and  dynamic  boundary  con- 
ditions at  the  surface  and  bottom  of  the  reservoir.  The  bottom  friction 
terms  are  represented  using  a Chezy  coefficient  in  the  following  form: 


where  C is  the  Chezy  coefficient.  The  terms  L and  L represent 

x y 

the  wind  shear  stress  effect  on  the  water  surface.  These  terms  are  of 


the  form 


and 


(7) 


T 

y 

(h  + n) 


where  T 

x 

surface. 


T are  the  wind  stress  components  acting  on  the  water 


10.  To  solve  the  governing  equations,  a finite  difference  approx- 
imation of  the  equations  and  an  alternating  direction  technique  are 


employed.  A space-staggered  scheme  is  used  in  which  velocities,  water- 


level  displacement,  bottom  displacement,  and  water  depth  are  described 
at  different  locations  within  a grid  cell  as  shown  in  Figure  2. 
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Figure  2.  Grid  system  and  variable 
definition  locations 

The  first  step  in  the  calculation  consists  of  computing  u and  n im- 
plicitly and  v explicitly,  advancing  from  time  nAt  to  (n  + l/2)At  . 
The  parameter  n is  an  integer  representing  the  time  step  at  which  the 
calculations  are  being  conducted.  The  second  step  computes  n and  v 
implicitly  and  u explicitly,  advancing  from  time  (n  + l/2)At  to 
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(n  + 1 ) At  . Central  differences  are  used  for  evaluating  all  derivatives 
in  the  governing  equations.  The  application  of  these  difference  approxi 
mations  gives  rise  to  corresponding  difference  equations  centered  about 
different  points  within  a grid  cell.  These  expressions  require  the 
evaluation  of  certain  quantities  at  locations  different  from  those  de- 
fined in  the  grid  system.  Such  quantities  are  replaced  by  values  com- 
puted from  a one-  and  two-dimensional  averaging  of  neighboring  values. 
The  time  interval  At  is  taken  as  the  time  required  to  complete  the 
full  cycle  in  the  computational  procedure;  however,  each  half  cycle  is 
treated  by  a different  set  of  equations  so  a system  of  six  operational 
equations  is  used.  The  difference  approximations  used  in  the  first 
step  of  the  alternating  direction  technique  are: 


n+1/2 


n+1/2  1 /At 


n+1/2 


n+1/2 


n+1/2  1 / At 


where  superscripts  refer  to  multiples  of  time  increments.  For  simplic 
ity  of  notation,  some  terms  have  been  maintained  in  differential  form 
within  angle  brackets  < > or  double  brackets  « ».  The  differential 
forms  are  defined  by  the  examples: 


11.  These  two  equations  can  be  solved  simultaneously  for  quanti- 
ties un+1/2  and  nn+1//2  along  a grid  line  k . The  additional  veloC' 
ity  component  vn+1/,£"  can  be  determined  explicitly  from  the  expression 


F*0 FET 


n+1/2 


At\  -n+1/2 


1 /At  \ n+1/2 


at  j,k+l/2 


The  equations  for  the  second  step  of  the  alternating  direction  tech- 
nique are: 


n+1  _ n+1/2 


,n+l  1 /At\-n+l/2 


n+1/2 


n+1/2  ^ 1 


n+1/2  at  j ,k+l/2  (13) 


n+1  _ n+1/2  1 /At 


n+1/2 


These  equations  are  solved  simultaneously  for  the  quantities  vn+l/2 
and  r|n+-'-/2  along  a grid  line  j . The  additional  velocity  component 
un+l/2  can  be  determined  explicitly  from  the  expression 


n+1/2  1 


n+1/2 


n+1/2  A 1 


n+1/2  at  j+l/2,k 


12.  The  formulation  of  the  problem  in  this  way  (ADI  method)  con- 
siderably simplifies  the  computational  procedure.  The  implicit  equa- 
tions are  solved  simultaneously  during  each  half  time  step.  However, 
with  this  formulation  scheme  the  coefficient  matrix  is  tridiagonal 

and  a relatively  simple  solution  procedure  is  available. 

13.  Two  types  of  boundaries  are  involved  in  the  calculations: 
the  solid  boundaries  at  fixed  coastlines  and  the  artificial  tidal  input 
boundaries  arising  from  the  need  to  truncate  the  region  of  computation 
(in  order  to  minimize  computational  time  requirements). 


jami  tma 
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F 

i 
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i . 


lU . A condition  of  complete  reflection  is  adopted  at  solid 
boundaries.  While  some  dissipation  does  occur  at  the  shoreline,  this 
should  not  be  significant  in  this  application.  The  actual  boundary 
condition  for  the  solid  boundary  can  be  written  as 

"V  = 0 (16) 

n 

where  V denotes  the  normal  component  of  velocity. 

15.  Artificial  tidal  boundaries  were  used  in  the  model  to 
describe  the  tidal  action  that  occurs  at  the  ocean  computational  bound- 
aries. These  boundaries  must  be  accurately  defined  since  the  tides 
applied  at  these  boundaries  represent  the  major  forcing  function  driving 
the  hydrodynamic  system.  The  water-surface  elevation  time  history  for 
the  desired  tidal  cycle  is  specified  at  each  such  boundary  and  applied 
during  the  operation  of  the  model. 

Application  to  San  Pedro  Bay 

16.  San  Pedro  Bay  is  formed  by  the  curvature  and  indentation  of 
the  southern  California  coastline  (Figure  3).  Sheltered  to  the  west  by 
Point  Fermin,  the  bay  is  open  to  the  south  and  southeast  except  for  the 
slight  protection  offered  by  Catalina  Island.  Originally  an  open  bay, 
the  protection  afforded  by  its  orientation  has  been  augmented  by  an 
8-mile-long*  breakwater  extending  from  Point  Fermin  eastward  to  near 
Seal  Beach. 

17*  The  breakwater  consists  of  three  sections.  The  San  Pedro 
breakwater  (oldest  of  the  three)  is  11,000  ft  long  and  extends  from  the 
shoreline  east  of  Point  Fermin  to  Angel's  Gate,  which  is  the  2,100-ft- 
wide  navigation  opening  for  Los  Angeles  Harbor.  The  Middle  breakwater 
is  18,500  ft  long  and  extends  from  Angel’s  Gate  to  Queen's  Gate,  which 
is  the  1,800- ft-wide  navigation  opening  for  Long  Beach  Harbor.  The 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Si)  units  is  presented  on  page  3. 
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Figure  3.  Site  map  and  limits  of  numerical  model 


Long  Beach  breakwater  is  the  third  section  of  the  breakwater  and  extends 
13,350  ft  due  east  of  Queen's  Gate. 

18.  A detailed  study  of  the  circulation  patterns  in  San  Pedro  Bay 
was  necessary  in  order  to  predict  the  effect  on  circulation  patterns  of 
changes  in  harbor  configuration.  To  get  the  required  resolution  of  the 
harbor  basins,  channels,  and  proposed  modifications  it  was  necessary  to 
use  a 300-ft  step  size  in  the  finite  difference  grid.  For  proper 
application  of  tidal  boundary  conditions,  the  region  of  study  was  ex- 
tended seaward  approximately  1.4  miles  from  the  Middle  breakwater. 

The  tidal  boundary  was  taken  as  parallel  to  the  Middle  breakwater  and 
the  region  of  application  of  the  model  is  shown  in  Figure  3.  A region 
of  about  60  square  miles  was  included  in  the  model  coverage  with  a fi- 
nite difference  grid  dimension  of  101  cells  by  199  cells  or  approximately 
20,000  finite  difference  cells.  There  were  approximately  12,000  water 
cells  where  the  velocity  components  and  tidal  elevations  were  calculated 
each  time  step  during  a tidal  cycle.  A time  step  of  45  sec  was  used  in 
the  numerical  model  so  a 25-hr  tidal  cycle  required  2000  calculation  time 
steps.  The  numerical  model  required  approximately  1-1/4  hr  of  CDC-7600 
computer  time  to  simulate  25  hr  of  prototype  time. 

19.  The  model  required  approximately  5 hr  of  prototype  time 
(15  min  of  computer  time)  to  "warm  up"  due  to  transients  associated 
with  starting  the  model.  Actual  calculations  for  a 25-hr  tidal  cycle 
required  approximately  1.5  hr  of  CDC-7600  computer  time  when  the  warm- 
up time  is  included . 

20.  Preparation  of  the  input  data  for  the  numerical  model  is 
extensive  since  each  of  the  20,000  finite  difference  cells  requires 
certain  input  data  concerning  the  cell  characteristics.  The  cell  must 
be  identified  as  a land  or  water  cell,  the  depth  of  the  water  relative 
to  some  datum  must  be  given,  and  the  friction  characteristics  of  the 
cell  must  be  identified.  The  friction  characteristics  are  expressed  in 
terms  of  Manning's  n which  vary  for  mud,  sand,  rock,  shoaling  areas, 
vegetation  in  marsh  areas,  etc. 

21.  After  the  input  data  have  been  prepared,  the  desired  tidal 
cycle  is  applied  to  the  model  and  the  results  are  examined  and 


compared  with  prototype  and  physical  model  results.  The  model  param- 
eters are  then  adjusted  until  the  desired  agreement  exists  between  the 
numerical  model  and  known  data. 


PART  III:  MODEL  VERIFICATION 

Existing  Conditions 

22.  Before  any  numerical  model  can  be  applied  as  a predictive 

tool  with  a reasonable  degree  of  confidence  it  must  be  verified.  It 

must  demonstrate  an  ability  to  produce  results  that  agree  with  known 

data  for  existing  conditions.  The  numerical  model  was  verified  by 

comparing  the  results  with  data  obtained  from  the  previously  verified 

2 

physical  model  of  Los  Angeles  and  Long  Beach  Harbors.  This  technique 
was  used  because  of  the  ready  availability  of  sufficient  data  of  the 
required  quality  and  quantity.  A typical  spring  tide  with  a 7.1-ft 
diurnal  range  was  one  of  the  tides  used  in  the  original  base  tests  on 
the  physical  model. 

Tidal  elevations 

23.  From  physical  model  tests,  surface  elevation  data  were 
available  at  13  locations  and  velocity  data  were  available  at  7 ranges 
(a  total  of  22  gages)  as  a function  of  time  for  an  entire  tidal  cycle. 
Locations  of  these  gages  are  shown  in  Figure  U.  In  addition,  photo- 
graphs of  surface  velocity  patterns  at  hourly  intervals  were  available. 

2h.  The  input  data  were  prepared  for  the  numerical  model  and 
the  7.1-ft  diurnal  spring  tide  was  applied  at  the  tidal  input  boundary. 
Results  were  compared  with  the  physical  model  data,  and  parameters  in 
the  numerical  model  were  adjusted  to  bring  the  numerical  results  into 
agreement  with  the  physical  model  results. 

25.  Tidal  elevations  calculated  from  the  numerical  model  were 
compared  with  physical  model  data  at  13  locations  distributed  through- 
out Los  Angeles  and  Long  Beach  Harbors,  and  results  are  shown  in 
Plates  1-5.  The  numerical  results  and  physical  model  results  are  es- 
sentially identical  at  all  gages.  This  result  coincides  with  past 
experience  which  has  indicated  that  tidal  elevations  are  reproduced 
extremely  well  by  the  two-dimensional  depth-averaged  numerical  model. 
Tidal  velocities 

26.  Reproducing  tidal  velocities  is  a much  more  difficult 
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ANGELES 


Figure  4.  Gage  locations  for  circulation  study 


problem  than  reproducing  surface  elevations.  This  difficulty  arises 
from  several  sources:  (a)  the  numerical  model  calculates  a depth- 

averaged  velocity,  an  entity  which  cannot  be  directly  measured  in  the 
prototype  or  physical  model;  (b)  prototype  and  physical  model  velocities 
cannot  be  measured  as  accurately  as  surface  elevations  due  to  instrument 
limitations;  and  (c)  the  numerical  model  calculates  velocities  at  dis- 
crete points  in  the  hydrodynamic  system,  therefore  the  finite  difference 
grid  locations  will  in  general  not  exactly  coincide  with  the  location  of 
velocity  gages  in  the  prototype  or  physical  model. 

27.  Despite  these  difficulties,  satisfactory  agreement  was  found 
to  exist  between  the  numerical  model  velocities  and  the  physical  model 
velocities.  Velocity  comparison  results  are  shown  in  Plates  6-27; 
numerical  model  velocities  are  shown  compared  with  the  velocities 
measured  at  the  midpoint  depth  in  the  physical  model.  The  agreement 

is  particularly  good  when  consideration  is  given  to  the  accuracy  of  the 
velocities  measured  in  the  physical  model.  The  low  current  speeds  ex- 
perienced in  San  Pedro  Bay  are  often  below  the  threshold  of  standard 
laboratory  meters  so  floats  were  used  to  measure  model  currents.  This 
measuring  technique  yields  a velocity  value  which,  rather  than  being 
an  instantaneous  value  at  a specific  point,  is  somewhat  averaged  over 
time  and  space.  Considering  the  possible  sources  of  error  in  the 
measurements,  the  physical  model  velocity  data  are  considered  accurate 
to  approximately  +0.2  fps. 

28.  No  attempt  has  been  made  to  draw  a smooth  curve  through  the 
numerical  results.  The  calculated  velocities  are  simply  plotted  each 
half  hour  for  the  tidal  cycle.  A smoothing  of  the  data  would  appear 
to  bring  the  numerical  results  into  even  better  agreement  with  the 
physical  model  data  since  the  numerical  model  picks  up  some  fluctua- 
tions in  the  flow  which  are  not  as  apparent  in  the  physical  model  data. 
This  is  particularly  true  in  the  Cerritos  Channel  gages.  These  rather 
short-period  fluctuations  in  velocity  would  tend  to  be  filtered  out  of 
the  physical  model  data. 

29.  The  overall  bay  circulation  patterns  were  obtained  from  the 
numerical  model  by  plotting  the  instantaneous  depth-averaged  velocity 


at  each  point  in  the  bay  as  a vector  at  specific  times  during  the  tidal 
cycle.  The  direction  and  length  of  the  vector  indicate  the  direction 
and  magnitude  of  the  velocity,  a vector  plot  is  illustrated  in  Figure  5. 
These  results  were  compared  with  surface  current  pattern  photographs 
from  the  physical  model  tests.  In  order  to  facilitate  a direct  compari- 
son of  results,  the  vector  plots  from  the  numerical  model  were  reproduced 

to  the  same  large  scale  used  in  the  surface  current  photographs.  The 
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vector  plots  have  therefore  been  published  as  a separate  report.  The 
overall  circulation  patterns  observed  in  the  physical  model  were  repro- 
duced very  satisfactorily  in  the  numerical  model.  This  is  particularly 
true  when  it  is  recognized  that  surface  currents  were  observed  in  the 
physical  model  while  the  numerical  model  produced  depth-averaged  veloc- 
ities. The  major  gyres  observed  in  the  physical  model  are  present  in 
the  numerical  model  results  and  extend  over  essentially  the  same  spatial 
region.  The  complete  extent  of  the  gyres  in  the  numerical  model  are  dif- 
ficult to  ascertain  when  the  vectors  are  plotted  to  the  scale  needed  to 
observe  the  major  circulation  patterns  in  the  entire  bay.  Thus,  the 
large  gyre  east  of  the  Angel's  Gate  opening  (Figure  5)  appears  to  be 
smaller  in  spatial  extent  than  that  observed  in  the  physical  model. 
However,  a larger  scale  plot  of  this  region  (Figure  6)  clearly  illus- 
trates the  actual  extent  of  the  gyre  to  be  very  similar  to  that  ob- 
served in  the  physical  model.  This  is  also  the  case  for  the  other  major 
gyres  although  in  some  cases  the  center  of  the  gyre  may  be  shifted  some- 
what from  that  observed  in  the  physical  model.  Any  apparent  differences 
in  the  strength  of  the  gyres  between  the  two  models  may  be  associated 
with  variation  in  the  strength  of  the  gyre  between  the  surface  strength 
and  an  average  strength  over  the  water  depth. 

Volumetric  flow  rates 

30.  The  volumetric  flow  rate  across  the  various  velocity  ranges 
was  also  calculated,  and  results  are  compared  with  physical  model  data 
in  Plates  28-32.  The  apparent  discharge  across  a range  was  calculated 
using  the  total  current  velocity  regardless  of  its  orientation  relative 
to  the  velocity  range;  the  adjusted  discharge  across  a range  was  calcu- 
lated using  the  velocity  component  normal  to  the  range  (Figure  7). 


Figure  5.  Velocity  vector  plot  of  tidal  circulation  patterns  for  existing  conditions 


VELOCITY  RANGE 


APPARENT  FLOW  ACROSS  AB  = (V)  (d)  ( h) 

ADJUSTED  FLOW  ACROSS  AB  = (VNORMAL)  ( d)  ( h) 

Figure  ?•  Apparent  and  adjusted  discharge 

31.  The  adjusted  discharge  will  obviously  always  be  less  than 
the  apparent  discharge,  the  ratio  of  the  two  quantities  being  the 
cosine  of  the  angle  0 . If  the  flow  across  the  velocity  range  is 
essentially  normal  to  the  range,  then  the  two  quantities  are  approxi- 
mately equal.  This  was  the  case  at  ranges  3,  5,  and  8 in  this 
investigation. 

32.  If  the  volumetric  flow  rates  are  integrated  over  the  entire 
tidal  cycle,  then  a net  flow  through  the  range  is  obtained.  The  net 
flows  obtained  in  the  numerical  model  and  from  the  physical  model  are 
compared  in  Table  1.  The  net  flows  are  small  numerically  in  comparison 
with  the  gross  volumetric  flow;  therefore,  small  percentage  errors  of 
variations  in  the  gross  flows  may  appear  large  compared  with  the  net 
flows.  For  this  reason,  the  net  flow  values  are  not  considered  exact 
and  should  be  used  only  to  indicate  flow  trends.  Net  adjusted  flows 
were  calculated  using  the  numerical  model  and  net  apparent  flows  were 
obtained  in  the  physical  model  which  is  primarily  responsible  for  the 
differences  observed  at  ranges  1 and  2.  The  primary  reason  for  the 
difference  observed  at  range  3 was  the  use  of  only  three  gages  to 
calculate  the  flow  across  the  range  in  the  physical  model  while  velocity 
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values  at  seven  locations  along  the  range  were  used  in  the  numerical 
calculations.  The  results  at  ranges  5 and  8 agree  extremely  well. 
Considering  the  indicated  qualifications  concerning  the  accuracy  of 
the  net  flow  calculations,  the  agreement  between  the  numerical  and 
physical  model  results  is  very  satisfactory. 

33.  From  the  numerical  model  there  appears  to  be  a net  westward 
flow  in  Cerritos  Channel  of  approximately  64  by  10  cu  ft/tidal  cycle. 
This  is  based  on  the  flow  across  range  8 which  should  be  more  accurate 
than  that  across  range  5 because  of  flow  conditions  crossing  the  veloc- 
ity range.  The  physical  model  results  indicate  a westward  flow  of 
approximately  50  by  10^  cu  ft/tidal  cycle  across  range  8.  There  is  a 
net  flow  through  the  breakwater  which  can  be  obtained  by  summing  the 
net  flows  across  ranges  1,  2,  and  3. 

Plan  1A-2 


34.  To  get  an  estimate  of  the  ability  of  the  verified  numerical 

model  to  predict  the  effect  on  harbor  circulation  of  modification  to 

the  harbor,  the  model  was  applied  to  the  master  plan  1A-2  configuration. 

The  only  changes  made  in  the  input  to  the  numerical  model  were  the 

changes  in  geometry  associated  with  plan  1A-2  as  shown  in  Figure  8. 

The  7-1-ft  diurnal  spring  tide  was  used  in  the  numerical  model.  The 

calculated  results  were  then  compared  with  physical  model  data  obtained 

3 

in  a previous  study. 

35-  Comparisons  of  tidal  elevations  are  shown  in  Plates  33-37. 

As  was  the  case  for  the  existing  conditions  comparison,  the  tidal  ele- 
vations are  essentially  identical. 

36.  Velocities  are  compared  in  Plates  38-55.  Basically  the  same 
type  velocity  agreement  was  obtained  for  plan  1A-2  as  was  obtained  in 
the  original  model  verification.  The  numerical  model  appeared  to  do 
a good  job  of  predicting  the  tidal  velocities  for  the  modified  harbor. 

37-  Vector  plots  of  the  overall  bay  circulation  patterns  for 
plan  1A-2  were  also  obtained;  a sample  plot  is  shown  in  Figure  9.  To 
make  the  changes  in  circulation  patterns  produced  by  the  modification 
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Figure  8.  Plan  1A-2 


more  apparent,  vector  plots  were  made  of  the  velocity  differences. 

These  plots  reflect  the  velocity  differences  in  the  bay  produced  by  the 

plan  1A-2  modifications  to  the  harbor;  a sample  plot  is  shown  in 

Figure  10.  Large-scale  vector  plots  for  plan  1A-2  and  for  the  velocity 

differences  produced  by  plan  1A-2  are  included  in  Report  2 of  this 
1* 

series  of  reports. 

38.  The  volumetric  flow  rates  across  the  velocity  ranges  also 
were  calculated;  results  are  shown  in  Plates  56-60.  When  the  plan  1A-2 
results  were  compared  with  results  for  existing  conditions,  the  follow- 
ing major  changes  were  noted: 

a.  Maximum  flow  rates  through  Angel's  Gate  (range  l)  were 
increased  by  about  10  percent. 

Is.  Maximum  flow  rates  through  Queen's  Gate  (range  2)  were 
decreased  by  approximately  50  percent. 

c_.  Little  change  was  produced  at  range  3 (east  end). 

d_.  Maximum  flow  rates  across  range  5 (Navy  Mole)  were 
decreased  by  about  25  percent. 

e^.  Maximum  flow  rates  across  range  8 (L.  A.  Channel)  were 
increased  by  about  20  percent. 

39*  Net  flows  across  the  velocity  ranges  were  calculated  and 
are  compared  with  the  physical  model  results  in  Table  2.  The  most 
significant  result  is  the  reversal  of  the  net  westward  flow  in  Cerritos 
Channel.  Both  the  physical  and  numerical  models  predicted  a small 
eastward  net  flow  in  Cerritos  Channel  for  plan  1A-2.  In  comparing  the 
net  flows  obtained  from  the  numerical  and  physical  models,  the  previ- 
ous statements  concerning  the  accuracy  of  these  numbers  must  be  con- 
sidered. These  net  flows  should  be  considered  only  to  indicate  flow 
trends. 

40.  Results  obtained  from  plan  1A-2  indicate  that  the  numerical 
model  is  capable  of  being  used  as  a predictive  tool.  It  will  apparently 
yield  essentially  the  same  results  as  those  which  would  be  obtained 
from  the  physical  model  and  can  now  be  used  to  investigate  the  proposed 
modifications  to  Pier  J. 


PART  IV:  ALTERNATE  PIER  J CONFIGURATIONS 

Modifications  Investigated 

Ul.  Two  basic  modifications  of  Pier  J,  STFP  2 and  STFP  3,  con- 
sidered for  a tanker  terminal  were  developed  after  preliminary  examina- 
tion in  a small  wave  tank  arid  are  shown  in  Figures  11  and  12.  Two 
cases  were  actually  considered  for  STFP  2,  a 300-  and  a 600-ft  opening 
in  the  eastern  end  of  the  breakwater  protecting  the  anchorage  area.  The 
three  modifications  indicated  as  STFP  2 (300),  STFP  2 (600),  and  STFP  3 
were  all  considered  for  a 7-1-ft  diurnal  spring  tide.  The  tide  was  ap- 
plied to  the  verified  model,  and  the  model  was  used  to  predict  circula- 
tion patterns  resulting  from  each  potential  modification  to  Pier  J. 

Results 


k2.  The  tidal  elevations  and  velocities  were  compared  with  those 
obtained  for  the  existing  conditions  at  the  various  gage  locations.  Ex- 
cept for  local  effects  around  the  Pier  J modification  configurations 
themselves,  STFP  2 (300),  STFP  2 (600),  and  STFP  3 all  produced  essen- 
tially the  same  overall  tidal  circulation  results  as  the  existing  condi- 
tion. Their  effects  upon  the  total  circulation  patterns  for  the  bay 
were  almost  identical.  The  tidal  elevations  and  velocities  produced  by 
STFP  2 (600)  are  presented  in  Plates  61-87  as  being  typical  of  all  three 
cases. 

h3.  Plates  6l-65  show  that  the  tidal  elevations  were  not  changed 
by  the  Pier  J modifications.  This  was  expected  since  the  modifications 
were  such  that  an  effect  on  tidal  elevations  should  result.  Plates 
66-87  show  that  the  velocities  at  the  gage  locations  were  changed  only 
slightly,  primarily  near  maximum  flood  and  maximum  ebb  conditions. 
Appreciable  changes  occurred  in  velocity  patterns,  however,  in  the  local 
Pier  J area  which  differ  somewhat  for  the  three  cases  investigated.  The 
vector  plots  of  the  circulation  patterns  for  STFP  2 (300),  STFP  2 (600), 
and  STFP  3 as  well  as  the  vector  plots  indicating  changes  in  circulation 
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Figure  11.  STFP  2 


patterns  from  the  existing  conditions  are  presented  in  Report  2 of  this 
series.4  Examples  of  these  vector  plots  are  shown  in  Figures  13  and  l4 
44.  The  volumetric  flow  rates  across  the  velocity  ranges  are 
presented  in  Plates  88-102.  Again,  all  three  cases  yield  essentially 
the  same  results  and  indicated  the  following: 

■a.  In  general,  only  a slight  difference  in  volumetric  flow 
rates  from  existing  conditions  except  near  maximum 
flood  and  maximum  ebb. 

b.  An  increase  of  approximately  12  percent  in  the  maximum 
flow  rate  through  Angel's  Gate  (range  l). 

£.  About  a 6 percent  decrease  in  the  maximum  flow  rate 
through  Queen's  Gate  (range  2). 

cl.  Very  little  change  in  the  flow  rate  across  range  3 
(east  end). 

_e.  A 10  percent  reduction  in  maximum  flow  rates  across 
range  5 (Navy  Mole). 

f_.  Very  little  change  in  maximum  flow  rates  produced  at 

range  8 (L.  A.  Channel)  although  there  appeared  to  be  a 
slight  change  in  phase  of  the  maximum  flow  rate. 


45.  Net  flow  rates  across  the  velocity  ranges  are  indicated  in 
Table  3.  Net  flows  across  the  velocity  ranges  did  not  change  apprecia- 
bly from  existing  conditions;  in  particular,  the  net  westward  flow  in 
Cerritos  Channel  still  appeared  to  be  retained  although  perhaps  reduced 
slightly. 


PART  V:  CONCLUSIONS 


L6.  The  Pier  J modifications  considered  in  this  investigation 
resulted  in  only  minor  overall  changes  to  tidal  circulation  in  the 
Los  Angeles  and  Long  Beach  Harbors  complex.  The  primary  circulation 
changes  are  local  in  nature.  All  three  modifications — STFP  2 (300), 

STFP  2 (600),  and  STFP  3 — produced  essentially  the  same  overall  effect 
on  tidal  circulation.  Specific  conclusions  concerning  the  changes 
in  tidal  circulation  produced  by  the  Pier  J modifications  are  as 
follows : 

a.  A slight  increase  in  flow  through  the  Angel's  Gate 
opening. 

b.  A slight  decrease  in  flow  through  the  Queen's  Gate 
opening. 

c_.  Very  little  effect  on  flow  in  Cerritos  Channel.  The  net 
westward  flow  in  Cerritos  Channel  appears  to  be  retained 
although  perhaps  decreased  very  slightly. 

cl.  Little  effect  northeast  of  the  Queen's  Gate  opening 
except  in  the  immediate  vicinity  of  the  actual 
modifications . 

e_.  Changes  in  the  surface  current  patterns  northwest  of  the 
Queen's  Gate  opening;  these  changes  are  primarily  confined 
to  the  Long  Beach  portion  of  the  harbor  complex  and  prin- 
cipally show  an  alteration  in  the  location  of  flow 
rather  than  volume. 
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Table  1 

Net  Discharge  per  Tidal  Cycle  for  Existing  Conditions 


net  flood  condition;  - = net  ebb  condition. 


Table  2 

Net  Discharge  per  Tidal  Cycle  for  Plan  1A-2 


Range 

Physical  Model 
Apparent  Net  Discharge 

10^  cu  ft 

Numerical  Model 
Adjusted  Net  Discharge 

10^  cu  ft 

1 

+1090 

+218 

2 

+80 

+698 

4 

3 

+1+80 

+223 

• 1 

5 

-l60 

+13.2 

8 

+80 

+31.1 

V ' 

Note: 


net  flood  condition;  - = net  ebb  condition. 


Table  3 


Numerical  Model  Adjusted  Net  Discharge  per 


Tidal  Cycle,  10  cu  ft 


Velocity  Range 
3 


Existing 

l+i+6 

61+1 

-35.3 

22.2 

-61+ 

Plan  1A-2 

218 

698 

223 

13.2 

31.1 

STFP  2 (300) 

1+70 

589 

131 

8.9 

-62.7 

STFP  2 (600) 

1+66 

61+3 

109 

22.1 

-61.9 

STFP  3 

1*45 

613 

152 

9.7 

-65.3 

Note:  - = net  ebb  condition. 
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APPENDIX  A:  NOTATION 
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Chezy  coefficient 
Coriolis  parameter 
Acceleration  due  to  gravity 
Water  depth 
Grid  lines 

Acceleration  effect  of  the  wind  stress  acting  on  the  water 
surface  in  the  x and  y directions 

Integer  defining  the  calculation  time  step 

Effect  of  bottom  roughness  in  x and  y directions 

Time 

Wind  stress  components  acting  on  the  water  surface 
u(x,y,z,t ) 

Depth-averaged  velocity  component  in  the  y direction 
Normal  component  of  velocity 

Depth-averaged  velocity  component  in  the  x direction 
Longitudinal  coordinate  measured  along  the  estuary  axis 
Rectangular  coordinate  variables 
Transverse  coordinate 
Vertical  coordinate 

Water-level  displacement  with  respect  to  datum  elevation 
p(x,y,z,t) 
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